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Abstract. 
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constraint operator in non-perturbative quantum gravity. In particular, we consider 
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exhibit the matrix elements of the hamiltonian constraint operator in the spin 
network basis in compact algebraic form. 
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1. Introduction 



In general relativity, the local dynamics of the gravitational field is governed by a 
constraint rather than by a hamiltonian Q . Accordingly, in non-perturbative quantum 
gravity the quantum dynamics is governed by a constraint operator Q , usually denoted 
the "hamiltonian constraint operator" (HCO) or the Wheeler-DeWitt operator. The 
construction and the analysis of this operator is a central problem in quantum gravity. 
For a long time, only formal or seriously ill-defined versions of the HCO were known. 
Concrete computations were thus impossible beyond semiclassical or minisuperspace 
approximations, and a predictive theory of the non-perturbative quantum dynamics 
of gravity was consequently lacking. 

The situation is perhaps improving with the development of the loop version 
of non-perturbative quantum gravity |3|], derived from the introduction of the Sen- 
Ashtekar variables j|, || (now evolved into the real Barbero variables [|[ 0). This 
approach, denoted "loop quantum gravity" , has provided a rigorous framework 
for describing quantum geometry (for recent reviews, and complete references see 
|||, |9[ [n]]). As far as kinematics is concerned, the approach has yielded striking 
predictions: in particular, geometric operators corresponding to volumes and areas 
have been shown to have computable discrete spectra. This result was first obtained by 
Rovelli and Smolin in | ll[] and later re-derived and sharpened with various alternative 
techniques ([^2[ O, |l4]; see |15| and references therein); its physical significance was 
first asserted in jig ]. 

As far as dynamics is concerned, an extensive array of investigations of the HCO, 
starting from the pioneering work of Jacobson and Smolin [ fL7[ and continuing in 
H [l^, |l9| 21 1, has recently culminated with the construction, due to Thiemann, 
of a mathematically well-defined operator [ p2[ . This construction is already being 
used in various directions of research, for instance in the search for exact solutions to 
the constraints J23[ , or as a key ingredient of a spacetime formulation of the theory 



j24, 25 . In this last context, the matrix elements of the HCO play a role analogous to 
the coefficients of the vertices in the Feymann diagrams. Modifications of the operator 
have been considered in [ pc[ . 

Thiemann's operator, however, is defined in a rather involved manner: its matrix 
elements and explicit action are far from easy to evaluate. This difficulty obstructs the 
research directions mentioned and renders the analysis of the operator quite difficult. 
In this work, we address this difficulty directly by studying Thiemann's operator in 
detail, and showing how its matrix elements can be computed explicitly. 

To perform the computation, we take advantage of two key techniques. First, 
we work in the spin network basis, introduced in quantum gravity in [^7[ pq| , 
and fully developed in |^9) (see |30| for an interesting historical reconstruction of 
these developments). Second, we use the technique developed in ||, based on the 
(Kauffman-Lins) tangle theoretical version of recoupling theory |3l[] (whose utility for 
quantum gravity was first suggested in |52]]). Since Thiemann's operator is gauge 
invariant but is defined in terms of gauge non-invariant operators, we extend here 
these techniques to gauge non-invariant states |53|, 

Furthermore, we take full advantage of the complete equivalence between the "loop 
representation" || and the "connection representation" |35|, [36| ^2| of loop quantum 
gravity (proven in detail in p7| ; see also [jl5| ||j|). It is this equivalence that allows 
us, for example, to use graphical calculus in a rigorous way in dealing with states that 
are functions of generalized connections. 
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We study here only the euclidean term of the HCO, but we expect that the 
techniques we use can be extended to the lorentzian term as well. Also, we restrict our 
attention to trivalent vertices: extension to any valence is straightforward. Thiemann 
defines two versions of the HCO: symmetric and non-symmetric. Here, we consider 
both operators, and show that they exhaust all "natural" orderings of the constraint. 

We do not discuss the physical correctness of Thiemann's operator. In particular, 
we do not address the objections that have been raised against it, concerning for 
instance the state dependence of the rcgularization, the strong commutativity (which 
shows consistency, but seems to fail to reproduce the classical constraint algebra [B8| ) 
or the possible need for adding terms to the operator suggested for instance in [p6|, |24| . 
Rather, we simply study the operator as defined by Thiemann, and show how its 
matrix elements can be computed, using recoupling theory. This analysis is likely to 
be needed before attempting a physical evaluation and, in any case, it is relevant even 
if modifications to the Thiemann operator are to be made. 

In section || we define the mathematical framework. In particular, we discuss the 
generalization of the spin networks to the gauge non-invariant case. In section ft| w e 
show in detail why graphical methods can be used in the present context, following |37| . 
Section |i| reviews Thiemann's construction of the hamiltonian constraint. Section 5 
presents our main calculations. For completeness, we have included three appendices: 
A discussion of the volume operator (Appendix A) ; technical details on the application 
of the graphical calculus in the context of the connection representation (Appendix 
B); and the basic formulae and concepts of recoupling theory (Appendix C). 



2. Loop quantum gravity. 

We start with a 4-manifold M. diffcomorphic to £ X R, where £ is a compact 3- 
manifold. The classical phase space of general relativity can be taken as the space 
of the canonical pairs (A l a ,Ef) of smooth fields over £. A\ is a connection 1-form 
on £ taking values in the Lie algebra of SU(2). Its conjugate momentum Ef is a 
vector density of weight one on E which takes values in the dual of su(2). We use the 
dimensional constants 

q _ 167rGNowton 



'0 — ll< - r — itm Plank — ^ i 

where GNewton is Newton's gravitational constant and Zpianck is the Planck's length. 
The fundamental Poisson bracket is: 

{Ai(x),E»(y)} = G6 i j 5 b a 6 3 (x,y). (2.1) 

The canonical variables are subject to three sets of constraints - the SU(2) Gauss 
constraint, the diffeomorphism constraint and the hamiltonian constraint. The space 
A of the SU(2) connections on E is the classical configuration space of general 
relativity @§. 

The quantum theory is defined as a linear representation of an algebra of classical 
observables on a Hilbert space 7i° . There exist now a number of constructions of 
the Hilbert space TL° needed for defining this representation (see Q), which have all 
turned out to be unitary equivalent. We will start here from the Schrodinger-likc 
representation of TL° as an L 2 space of functionals over the configuration space Q . 
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More precisely, as it is usually the case in quantum field theory p9[ , TL° is defined 
as a space of functionals over a suitable completion A of the configuration space A, 
formed by "generalized connections" , or "distributional connections" . 

2.1. State space 

In order to define a generalized connection, we need some definitions. An edge e is an 
oriented, 1-dimensional sub-manifold of E with two boundary points, called vertices, 
which is analytic everywhere, including the vertices. A generalized connection A in A 
is defined |l3|, [hj) as a map that assigns an element A(e) of SU(2) to each oriented 
edge e in £, satisfying the following requirements: i) A(e^ 1 ) = (^4(e)) _1 ; and, ii) 
A(e2 o ei) = A(e2) ■ A{e\). Here e _1 is obtained from e by reversing its orientation, 
&2 oei denotes the composition of the two edges (obtained by connecting the end of ei 
with the beginning of and A(e2) • A(e\) is the composition in SU(2). Any smooth 
connection A is also a generalized connection: the group element that it assigns to 
an edge e is the parallel transport along e, that is A(e) = he(A) := Pexp(— J„ A). 
Accordingly, we refer to A(e) as "the parallel transport of the generalized connection 
A", in the same spirit in which we write 6(f) = J dxS(x)f(x) for the Dirac 6- 
"function" . One can view the space A of the generalized connections as the closure of 
the configuration space A in a suitable topology. 

There exist a natural measure fx on A, induced by the Haar measure on SU(2), 
denoted the Ashtekar-Lewandowski-Baez measure. Its explicit form can be given as 
follows. A graph 7 in £ is a collection of edges such that if two distinct edges meet, 
they do so only at the vertices. A vertex v is called n-valent or of valence n if it has 
n adjacent edges. Fix a graph 7 with N edges ex, ...,ejv, and consider the functions 
\P<y,^ on A of the form ^ 1 ^(A) = ^(A(ei), A(ejy)) for some smooth function tp 
on [SU(2)] N . These are cylindrical functions, since they depend only on a finite 
subset of the infinite "coordinates" of A. The measure \i° is defined by (its cylindrical 
projection) 

[_d fJ ,°(A)%^(A)= / d n f x H (9i,...,9N)^(gx,...,gn), (2.2) 

J A JsU{2) n 

where [in is the Haar measure on SU(2). We use the Hilbert space H° :— L 2 (A,d/i°) 
of square-integrable functions on A as the (kinematical) Hilbert space of the theory. 
( "Kinematical" in the sense of Dirac - the physical Hilbert space is then obtained 
by imposing the quantum constraints.) The elements of TL° can be viewed as 
wave functions of (generalized) connections, analogous to the wave functions of 
distributional fields in the Schrodinger representation of a free field theory. The 
cylindrical functions are dense in 7i°; they represent "typical" normalizable states. 

Let us now discuss the SU(2) gauge invariance. A generalized gauge transformation 
is a map g which assigns an SU{2) element g(x) to each point x of E. It acts on A 
at the end points of edges: A(e) — * t;(t; + )~ 1 • A(e) ■ g(v-), where w_ and v + are the 
beginning and the end point of e. The space of the equivalence classes of generalized 
connections under generalized gauge transformations is denoted A/Q. The measure 
li° induces a natural measure ft° on A/Q , obtained by the push-forward of fj,° under 
the projection map that sends A to A/Q. 

The gauge invariant functions form the SU (2) gauge invariant Hilbert space 
H" := L 2 (A/Q ,djl ). These are the solutions of the SU(2) constraint. Since the 
spaces under consideration are compact and measure-normalized, we can regard H° 
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as the gauge invariant sub-space of the Hilbert space TL° := L 2 (A,dp°) |l3|, [hj (4l[ . 
Thus, the SU (2) gauge invariant quantum states can be expressed as complex- valued, 
square-integrable functions on A/G, or, equivalently, as C?-invariant square-integrable 
functions on A. 

2.2. Spin network basis 

In order to work on a Hilbert space it is convenient to have a basis. A very convenient 
basis in H.° is formed by the spin network states. These are obtained by using 
the fact that functions on SU(2) can be expanded in irreducible representations. 
Given N irreducible representations ji, ■■■jjjsr of SU(2), we define an invariant tensor 
c mk+1 — mN mi ... mk as a multi-linear map from ® J=1 ij to ^>f =k+1 ji that transforms 
covariantly, namely such that: 

c nk+1 - nN ni ...n k = (2.3) 

for arbitrary g S SU(2), where j/(<7) is the matrix representing g in the representation 
ji. Such an invariant tensor is also called an intertwining tensor from the 
representations ji , . . . , jk to , . . . , ■ All invariant tensors are given by the standard 
Clebsch-Gordan theory. 

An extended spin network s is defined as a quintuple s = (7, j, c, /?, M ) consisting 

of: 

• A graph 7 in S; 

• A labeling j :— (ji, .., j N ) of the edges ei, ...,e N of the graph 7 with non-trivial 
irreducible representations of SU(2). We refer to pi, equal to twice the spin of 
the representation jj, as the color of the edge ej; 

• A labeling p :— (pi, py) of the vertices V\,...,Vv of 7 with (possibly trivial) 
irreducible representations of SU(2), with the constraint that for every vertex v a 
the representation p a sits in the tensor product of the representations assigned 
to the edges adjacent to v a ] 

• A labeling c = (ci, ...,cy) of the vertices v±, ...,vy of 7 with invariant tensors; 
more precisely, assigned to a vertex v a is an intertwining tensor c Q from the 
representations assigned to incoming edges and p a to the representations assigned 
to the outgoing edges; 

• A labeling M := (Mi, My) of the vertices Vi, ...,vy of 7 which assigns a vector 
M a in the representation p a to every vertex v a . 

The (extended) spin network state ^> s is defined as follows. Consider the "parallel 
propagator" A(ei) of A along each edge e/ of 7, in the representation jj. Contract 
these parallel propagators at each vertex v a , using the invariant tensor c a . Contract 
the leftover indices of c a at each vertex v a with the vector M a assigned to the vertex. 
Formally: 

N V 

*s(A) := [(g) j/(i(ej))8gM a ] • [®l =1 c a ]. (2.4) 

1=1 a=l 

Here '•' stands for contracting, at each vertex v a of 7, the upper- indices of the matrices 
corresponding to all the incoming edges, the lower indices of the matrices assigned to 
all the outgoing edges and the upper index of the vector M a with all the corresponding 
indices of c a . "i> s is a C°° cylindrical function on A. 
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The introduction of the vectors M has the sole purpose of avoiding the choice of 
a preferred basis in the SU (2) representations. If we introduce such a a preferred 
basis, we can replace the M's with the specification of a basis vector. In other words, 
the extended spin network states have a free index at each vertex v a , living in the 
representation p a . Below, we will find it more convenient to use the physicist-style 
indices, rather than the mathematician-style M notation. 

For any spin network s, ^ s is a function on A which is square-integrable with 
respect to the measure pf. We can fix a basis in the (finite dimensional) space of the 
invariant tensors at each vertex, and a basis in each representation p a , and restrict c 
and M to basis elements. The resulting set of spin network states form a basis in 7i° . 

It is easy to see that a spin network state ^ s is gauge invariant if and only if 
all the representations in p are trivial. The gauge invariant spin network states are 
therefore characterized solely by the graph 7, the labelings j of the edges and the 
intertwiners c at the vertices. The inclusion of the representations p and the vectors 
M at the vertices yields the extension of the spin network technology to the gauge 
non-invariant states, namely the extension from H.° to 7i°. 

2.3. Loop states 

An open loop I = (ei, . . . e„, M_, M + ) is defined by an oriented open line in E formed 
by a sequence of edges ei, . . . e„, and by two vectors (M_, M+) in the fundamental 
representation of SU(2). A loop a = (e 1: . . . e„) is defined as a closed line formed 
by a sequence of edges. An open loop state is defined by ^i(A) = M_ ■ j(A(e 1 j) ■ 
... ■ j(A(e n )) ■ M+. A loop state is defined as *f> a (A) = Tr\j(A(e 1 )) ■ ... ■ j(A(e n ))], 
where j is the fundamental (spin 1/2) representation of SU(2). Sums of products of 
open loop states span Tt°, and sums of products of loop states span H°. Products 
of loop states are called multiple loop states. Thus the multiple loop states span fi°; 
however, they do not form a true basis, but only an overcomplete basis. This is the 
overcomplete basis that was used in loop quantum gravity before the introduction of 
the spin network true basis. 

The relation between spin network states and loop states is as follows. Standard 
SU(2) representation theory shows that the spin network states can all be obtained as 
linear combinations of products of open loop states. This is because all representations 
of SU(2) can be obtained by tensoring the fundamental representation. Explicitly, 
these linear combinations are obtained by "symmetrizing the loops" along each edge 
as illustrated in P7|| . Conversely, every product of fundamental representations can 
be decomposed in irreducible representations, and therefore multiple loop states can 
be expanded in the spin network basis. Thus, the decomposition of quantum states 
in the spin network basis is the decomposition of functions on SU (2) in irreducible 
representations. See |9[ |7) for details. 

2. 4- Operators 

We shall not give here a complete construction of the operator algebra defined on 7i°, 
but only a few remarks on the operators that enter the Thiemann's definition of the 
HCO. These are the parallel transport operator and the volume operator. For each 
edge e, the parallel transport operator h[e] (with a free spin 1/2 index at each end 
of e) represents the classical quantity h e (A) := V exp(— J A), that is, the parallel 
propagator of the connection along e. This operator is diagonal, and is given by 
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multiplication by A(e): 

h[e] <3>j (A) = A(e) ^f(A). (2.5) 

The action of this operator on the spin network basis is not immediate. We obtain 
it by decomposing th e st ate /ifej'l's in the spin network basis, and this has to be worked 



out explicitly using (2.4). For instance, if e and the graph 7 do not intersect, then we 
have immediately that ft,[e]^ / s = $( sUe ), where (s U e) is the spin network obtained 
by adding the (disconnected) edge e to 7, and coloring it and its two open vertices 
with the spin- 1/2 representation. As a second example, assume that e sits entirely on 
the interior of one of the edges of s, say with label j. The decomposition of ft[e]W s 
in spin network states can then be obtained using the fact that the tensor product of 
two SU(2) representations is a sum of irreducible representations. Since A(e) is in the 
1/2 representation, we have in this case 

h[e]^ s = a+* s+ + a_# s _ (2.6) 

where s± are the spin networks obtained in the following way. Add two vertices to s, 
in correspondence to the end points of e. Label the edge between these two vertices 
by j ± 1/2. The representation p on the two new vertices is 1/2, and the intertwiner 
is the unique one on the three representations j,j ± 1, 1/2. The free indices on the 
new vertices are the free indices of the operator, and, finally, the coefficients a± are 
the two Clebsch-Gordan coefficients of the expansion of j <E> 1/2. (See appendix C.) 
More complicated cases can be worked out in a similar way. 

The rest of the elementary operators defined on Ti° contain Ef. These can be 
defined using the T n loop operators technology || , or using left-invariant vector fields 
on the group [ |35| , [36|. Thiemann's HCO depends on the volume, which is a composite 
operator based on such elementary operators. 

The volume operator has been constructed by Rovelli and Smolin in JlTf using 
the T n loop operators technology. This construction is based on a regularization and 
quantization of the classical volume. (An important clarification role was played by 
Loll's work [|l2| using a lattice regularization p2|.) Later, Ashtekar and Lewandowski 
in |l| used the left-invariant vector fields technology to construct a "mathematically 
natural" operator Val ■ They presented it as "closely related" [jl3] to the construction 
of Rovelli and Smolin and to the classical volume. Lewandowski has clarified the 
relation between the two operators, showing that the difference between the two can 
be interpreted as due to different regularization choices - not to the technique employed 
(both operators can be constructed in each approach) jl5|. See Appendix A. The two 
operators have the same action on generic trivalent vertices (tangents of adjacent edges 
non coplanar), which is the only case we consider here. In the following we need the 
explicit action of the volume operator. This has been worked out by DePietri and 
Rovelli in [^). See Appendix A. Thiemann has also derived various matrix elements of 
the volume operator in [43] using a different technique and confirming the results of 



3. Graphical binor representation 

In |37[ it was shown that computations in the connection representation of quantum 
gravity can be performed using the graphical Penrose binor calculus fl44| (as originally 
suggested by Smolin). The main ingredients of this calculus are summarized in 



(A) 2 
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Figure 1. A spin network states and its graphical representations. 



Appendix B . Furthermore, it was shown in |p7| that this graphical binor calculus is 
fully equivalent to the loop representation developed in || [ll|, |27], |j| . These results were 
obtained in the gauge-invariant context, but there is no restriction to gauge invariant 
functions in the rules of Penrose binor calculus, therefore they extend immediately to 
the case of gauge non-invariant cylindrical functions. Let us review these results. 

Following [[37], we associate a graphical representation a' la Penrose to each 
cylindrical function vP^f. As a first step, we considers a projection of the graph 7 
on a plane and construct its extended planar graph T ex (see figure |l|(B) and The 
extended planar graph is a two dimensional "thickening out" the graph: an oriented 
two-dimensional surface with the topology of 7X [0, 1]. Next, we expand the cylindrical 
function as a sum of products of irreducible representation on the edges. For each term 
of the sum, we insert the Penrose graphical representation of the matrices ji(A(ei)) in 
the corresponding edges er of the extended planar graph, using equation (B.6). Next, 
we insert the graphical representation of the contractors c Q in each vertex v a , using 
equations (B.7) and ( |B.cf ). Finally, lines corresponding to summed dummy indices 
are connected, with an additional factor (— l) p for each trace (closing a loop); here 
p is the color of the representation in which the index lives (see ( |B.3[ )). It follows 
from equation (B.9) that the planar graphical representation is independent from the 
orientation of the edges. Finally, we notice that the Penrose convention of indicating 
tensors by means of filled boxes is superfluous in this context, because a colored line 
that runs along an edge of the extended graphs always carries a parallel propagator. 
Thus the information given by the boxes in already contained in the fact that the line 
runs along an edge, and we can drop the boxes. 

In particular, we obtain in this manner a graphical representation for each spin- 
network state. Let us focus on these. First of all, in a spin network state we have 
a single irreducible tensor on each edge. Therefore each edge carries a single colored 
line (recall: color is twice spin). Second, at each vertex v a we have an intertwining 
tensor c a . The graphical representation of each invariant tensor can be given by means 
of a trivalent tree-like graph (with open ends connected to the lines in the adjacent 
edges). This is a consequence of equation flB.Sj ). We denote the "internal" edges of 
this trivalent decomposition as "virtual edges" . It is important to keep in mind that 
virtual edges are drawn with finite length on the extended planar graph, but they do 
not represent finitely extended objects on the space manifold S - they are simply a 
graphical representation of the index pattern of combinations of the SU{2) indices of 
Clebsch-Gordan coefficients. They characterize an intertwiner by giving its expansion 
in Clebsch-Gordan coefficients. Finally, at each vertex we have a representation p a , 



f More precisely we define the graphical binor representation of polynomial cylindrical functions, i.e, 
functions that are polynomial in the SU(2) matrices ge = he (A) associated to each edge e of 7. 
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and a vector M a (or a free index) sitting in this representation, contracted with the 
invariant tensor. Accordingly, in the graphical representation of the spin network 
state we have one of the virtual edges which is not attached to any external line. This 
free end is missing if p a is the trivial representation, namely if the vertex is gauge 
invariant. 

In figure [l] we show, as an example, (A): a spin network (in the connection 
representation) determined by an oriented graph, and the labels of edges and vertices. 
(B) A spin- network in the loop representation and its extended planar graph (as 
in |^). (C) The planar binor representation on the extended planar graph of the 
cylindrical function given by the spin-network. 

A simple advantage of this graphical representation is due to the fact that planar 
graphs are easier to manipulate than expressions with a large number of indices. 
But the main advantage is due to the peculiar nature of the tensors involved, which 
are based on SU(2) matrices. These tensors satisfy some basic relations that are 
translated into simple graphical formulas. The prototype of these is equation (B.4), 
from which (as shown by Kauffman and Lins pH ) all the other graphical relations can 
be derived. These graphical relations are the reason for the great utility of graphical 
recoupling theory, well known and much exploited in nuclear, atomic and molecular 
physics (where 50(3) plays a major role). The important formulas satisfied by the 
graphical representations of the quantum states are summarized in Appendix C. The 
calculations in this paper rely heavily on these formulas. 

As an example, consider the action of the holonomy operator h[ei] on a spin 



network state i& s , given in (2.5). Consider the particular case in which e sits on 
the interior of an edge e' of s. Following the prescription given above, the graphical 
representation of the r.h.s. of (|2.5|) is simply obtained by adding a line colored 1, 



inside the (expanded) edge e'. The expansion of this state in spin network states is 



then immediately given graphically by (C.1C) 



Notice that, following |37[, we have fully recovered in this way the "pure" loop 
representation of quantum gravity, in which states and operators can be defined 
without any reference to the Schrodinger-like functions of connections. In the same 
fashion, the Dirac energy-basis \n) representation of the harmonic oscillator allows us 
to formulate and solve the theory without reference to its Schrodinger representation. 
However, here we have recovered the loop representation starting from the gauge 
non-invariant connection representation. Thus we have obtained the generalization of 
the loop representation to gauge non-invariant states. This generalization is simply 
contained in the fact that the spin networks can have (colored) open ends at the 
vertices. Of course, one might have avoided the detour through the connection 
representation, and might have directly defined the loop representation as a suitable 
representation of an appropriate algebra of gauge non-invariant observables. This can 
be done, but we have preferred here the long way, in order to make the link with the 
connection representation (using which Thiemann constructed the HCO) fully explicit. 



4. Thiemann's Hamiltonian constraint operator. 

The Hamiltonian constraint of general relativity can be written as 



Vdet(g) 
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Figure 2. An elementary tetrahedron and the orientation choice of the edges. 



TL 



Vdet(g) 



Tr[K a ,K b ][E a ,E b 



(4.1) 



F ab is the curvature of A\ L , and K a is related to the extrinsic curvature of the spatial 
manifold S. This is the form that the constraint of the Lorentzian theory takes when 
written in terms of the real connection. Both (A^^Ef) and (Ef,Kl) are canonical 
pairs of variables. K\ is related to A\ by A\ — T l a — sgn(e)iQ, where T l a is the 
spin-connection annihilating the triad and sgn(e) is the sign of the determinant of 
undensitized triad e % a . There are two important differences between this connection 
and the one introduced by Ashtekar. First, all variables are real, while the Ashtekar 
connection, A\ — T l a + i sgn(e)-£Q, is complex. Second, the constraint (4.1) contains 
the factor of one over the square root of the metric. In Ashtekar's theory this factor is 
absorbed into the lapse function. The first term in the right hand side of (4.1) defines 
the hamiltonian constraint for gravity with euclidean signature - accordingly, 7i E is 
called the euclidean hamiltonian constraint. 

In order to simplify the hamiltonian constraint, Thiemann [22]] uses the following 
identities 



[E a ,E 



b] i 



^/detq 



G 



e abc {A l c ,V} 



and 



(4.2) 



G 



{Ai,K}, 



where V is the volume of X: 



V = j d 3 xy/detq = j d 3 x^ 



e abc eV k E?E*El 



g| -iwc' w j j k 

and K is the integrated trace of the densitized extrinsic curvature of S: 



K = I drx Kl E? 



(4.3) 



(4.4) 



(4.5) 



Using the identities (4.2) and (4.3) we can write the euclidean hamiltonian constraint 



as: 



U E [N] = ^Jj Z * t abC Tr(F ab {A c: V}), 



(4.6) 
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and the total hamiltonian constraint (4.1) then takes the form: 

H[N] = H E [N] -^J^ d 3 xN(x)e abc Tr ({A a , K}{A bl K}{A C1 V}) .(4.7) 



Another useful identity expresses the extrinsic curvature (4.5) in the form: 

K = ±{V,H E {1}}, (4.8) 



where 7i E [1] is the euclidean hamiltonian constraint from ( |4.6|) for lapse function set 
to one. 

By plugging ( |4.6| ) into (4.8) and then the result into ( |4.7|) we have an expression for 



the lorcntzian hamiltonian constraint defined entirely in terms of the volume and the 
connection. It is then relatively simple to turn this expression into an operator. In the 
connection representation, the connection A l a and the curvature F* b act (essentially) 
multiplicatively. The differential part of the HCO, originating from the action of 
the triads, is captured by the volume operator. But the volume operator has a well 
understood action on the kinematical state space. 

From now on, we concentrate on the euclidean term of the constraint. In the 
quantum theory we do not have an operator correponding to the connection, but only 
an operator corresponding to its parallel transport. Accordingly, we regulate Ti. E , 
expressing it in terms of parallel transports. We introduce a triangulation T of S 
formed by elementary tetrahedra A with analytic edges. We denote by sf(A), or 
simply by Si,i = 1,2,3 the three edges meeting at a vertex v(A) of a tetrahedron A. 
Also ay (A) := s;(A) o s y (A) o sJ 1 (A) is the loop, based at v(A) and composed out 
of the edges s,-, sj 1 , and the corresponding edge of the tetrahedron (see figure ||). 
To start with we write the Hamiltonian constraint as a sum over all tetrahedra in the 
triangulation 

H E [N] = J2 nl[N] = J2 I d " x N ( x ) £QbCTr (F ab {A c , V}) . (4.9) 
Aer Aer^ A 

In the limit in which a tetrahedron A shrinks to a point v(A) the Hamiltonian 
constraint over A is 

n%[N] = N(v(A)yi k Tr(h atj h Sk {h- k \V})+0(V(Ay/ 3 ) (4.10) 

where V(A) is the volume of the tetrahedron A. In the above expression h ai . is the 
holonomy along the closed loop ctij and h Sk is the holonomy along the open loop 
segment s^. That the above limit gives the correct expression (^1]) follows from the 
expansions: 

lim ft =l + lF ab sts b j (4.11) 

and 

lim h Sk = l + A c s c k . (4.12) 

A—>v(A) 

Thus, if we have a sequence n <— > T n of finer and finer triangulations of E, we can 
write 

H E [N] = lim V Ha[N]. (4.13) 
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Figure 3. The "small open loops" of the Hamiltonian operator. Notice that in 
gauge non-invariant states (states with open ends in the graphical representation) 
orientation matter. Two open ends can be joined only if their orientation is consistent. 



Since the operators corresponding to the classical quantities h and V are 
well understood, this equation provides an expression for the constraint which is 
appropriate for quantization. However, by simply turning this expression into an 
operator, with an arbitrary chosen sequence of triangulations, one obtains meaningless 
results. Following Thiemann, we define the quantum Hamiltonian operator by 
adapting the regularization to the state. For every state ^j t ip, we define 

H E [N] = lim Y Hi[N] * 7 ,^ (4.14) 
Aeri T) 

where 

Hl[N] = N(v{£))Hl = N v Hl 

Hi = - J^Tr J''' lu hJ.V^j , (4.15) 

and T„ (7) is a sequence of triangulations satisfying the following condition: (for each n) 
there is a vertex i>(A) sitting on each of the vertices v of 7; and if the vertex v is fc-valent 
then the corresponding vertex w(A) has precisely k adjacent edges, each contained 
in one of the k edges of 7 adjacent to v. The triangulation is arbitrary elsewhere. 
Since the classical expression converges to the correct hamiltonian constraint for any 
sequence of triangulations, this choice does not spoil the formal classical limit of 
the operator. On the other hand, this choice yields a well defined, diffeomorphism- 
covariant and non-trivial operator. This completes the definition of the quantum 
operator. 

Simple manipulations then show that the action of this operator can be written 

as 

H*[N]% t i,= 8N - E # a^v*7,^ (4-16) 

ueV(7) v(A)=v ^ ' 

where V(s) is the set of the vertices of s; is defined by (4.15); pa is one, whenever 
A is a tetrahedron with three edges coinciding with three edges of the spin network 
state joined at the vertex v, and zero otherwise. E(v) is the number of non-coplanar 
triples of edges of 7 at the vertex v. In the case of an n-valent vertex having no triples 
of coplanar edges E(v) = n(n — l)(n — 2)/ 6. Also N v is the value of the lapse function 
at the vertex (see Q for details). 
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Clearly of (4.15) it is not the only operator that has (4.10) for classical limit. 
Any other ordering of the operators h mj , h Sk and [/i" 1 ,^] gives ( 4.10D . We will 
consider the action of Thiemann's Hamiltonian constraint with the two orderings 



n 



Ea> 



— e ijfe Tr 
3iL 2 



n 



E(2) 



hs h [K k \V] 



(4.17) 
(4.18) 



which correspond to the cases analyzed by Thiemann j2i| . 

To compute the action of (4.17) and (4.18) on trivalent vertices, it is sufficient to 
compute 



Ha® *7,V> 



2i 






e ljkr Tr 


3lf 




2i 






e ljk Tr 





h Sk Vh Sk 



*7,V> 



(4.19) 
(4.20) 



because the volume annihilates the gauge-invariant trivalent vertices and therefore 
the other term of the commutators vanishes. In section 5.3, we show that these two 
choices exhaust all possible "natural" orderings. 



5. Calculation of the HCO matrix elements 

The action of the HCO is local, namely it is a sum of independent actions on each 
vertex. Therefore we can focus on its action on a single vertex v. Let us consider the 
action of the operator fi^ on a single trivalent vertex with edges e i; e^, e^ with colors 
p, q, r. We indicate the part of the spin network state, containing only the vertex v 
and its adjacent edges as | v). 



5.1. First ordering choice 

We investigate the action of the Euclidean Hamiltonian operator 7i^ a) on | v) 



H E A m I v) 



2i 

'3? 



ijk 



h[c 



h[c 



.in 



h[s k ]Vh[s^] | v) 



(5.1) 



As described in sections 2 and 3, /j[Sfe ] attaches an open loop segment s^T 1 to the 
edge efc with color r. It creates one additional vertex one^ and alters the color of 
the edge segment between the two vertices of Sfc. From ( C.lOp we have 




e=±l 




(5.2) 



From Appendix C, the coefficients a±(r) are 



a+(r) = 1, a-(r) = - 



r + 1 



(5.3) 
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Notice that the (lower) open end in the right hand side of (5.2) is located at the 
vertex. Thus, the vertex has now become gauge non- invariant: it carries the spin-1/2 
representation, with a free index, and its interwiner is determined by the the virtual 
edge of color r. 

The volume operator acts on this vertex. The action of the volume operator has 
been computed for gauge invariant states only, but a moment of reflection shows that 
the results in the gauge invariant case are relevant for the present case as well. In 
fact (in both its definitions) the volume operator acts on the (real!) edges adjacent 
to the vertex, giving a sum of terms - one term per triple of edges. If we act on a 
gauge-invariant 4-valent vertex, we obtain a sum of four terms, one for each triple of 
edges. In the present case, the volume acts on the three real edges, giving a single 
term, but the action of this single term is precisely equivalent to the action of one of 
the four terms of the volume over a four-valent gauge invariant vertex (See Appendix 
A). The free index of the vertex behaves as the fourth edge, but it is "invisible" for the 
operator. Therefore the action of V can be immediately obtained from the expression 
of V on 4-valent gauge- invariant vertices, which was given in ||, and is recalled in 
Appendix A. In particular, we show in Appendix A that for the combination of colors 
in (5.2) this action is diagonal. Using (A. 11), we obtain 




where w( , , , ) is given in (A. 10). 

The next step is to compute the action of /i[aiy]/i[sfe] on the right hand side of 
(5.4). This operator attaches ay and to the open ends in (5.4) according to the 
prescription given by the orientation of the open loops. Finally, the trace connects the 
free indices. This is represented in the graphical computation by tying up the open 
ends of the virtual edges (see figure f|). Moreover, we have to add a minus sign due 
to the trace prescription of the graphical representation (see section ||). We obtain 



~h[aij]h[s k ] - h[aji]h[s k 
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£=±1 e=±l 




(5.5) 



In the last equation we used the fact that the anti-symmetrization of tw o loops is 
equivalent to two 3-vertices of type (1,1,2) and we applied equation ( C.1C| ) to the p 
and q edges. Now, the tangle around the vertex in the right hand side of (5.5) can be 
evaluated using recoupling theory. We have (see Appendix C) 





® 




(5.6) 



V + e q + e 



In (:k6) we have not shown the bottom link of (5.5). Using equations (|C.ll ) end (C.9) 
we expand the right hand side of (5.6) 

T I 



p + e q + e 





(5.7) 



p + e q + e 



The chromatic evaluation of the closed nets in equations (|5.6|) and (5.7) gives 
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/ p-q 
2r 

(q-p + r)(2+p + q) 

4r(l + g) 
(p-q + r){2+p + q) 

4r(l+p) 
(p - q){p + q - r){2 + p + q + r) 



= < 



8r(l+p)(l + g) 



e = +l,e = +1 
e = -l,e = +l 
e = +l,e = -1 
e = — 1, e = — 1 



We have obtained 



Tr 



where 



ft™,. 




Z o I] A (1) (p,e;g,?;r) p+E 



e,e=±l 




-4 (1) (P, e;g,e;0 = X] I V^M^+M) #i(r,e) #2(p, e; g, e, r). 



(5.10) 



e=±l 



Collecting our results we obtain that the action of Thiemann's Euclidean HCO on 
a generic 3-vertex \v(p,q,r)} is given by: 



*2 B Kp,?,r)> = ^5: £ 



e=±l g=±l 



^ (1> (p, e; g, e;r) 



P 



+A (1) (g,e;r, e;p) 




+ A (1) (r, e;p, e; g) 




(5.11) 



where the coeficients A 1 are given by (5.10). 



5.2. Second ordering choice 

In the case of the ordering choice given by (4.20) we have to perform the computations 
in the reverse order. First we compute the action of the elementary loop operator 
ft" 1 [h aij — h aji ] on | v), and then, the action of h Sk V. Using the same technique as 
in the previous section we get: 
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= E E E a e{r)a E {p) ai {q) 



e=±l e=±l e=±l 



E X! a e (O a e(p)ag(g) j ^ 

e,e=±l J=r,r+2e ^ 




E E «-(r){; 



r+E 



r r + e J 
2 1 



jk(p, e;g,e; r, J) 



P+E 



•*" 1 



We first inserted the loop corresponding to the h[...] operator in the extended 
planar graph corresponding to the 3-vertex. Then we decomposed the result, using 
equation ( p. 10 ). Then we have applied the recoupling theorem ( p.6[ ) to the 4- vertex 
(r, r + e, 1, 2). Finally we have reduced the 3-vertex (p + e, q + e, J) to its fundamental 
form as in (5.5). In the last step we have defined: 




R2(p,e;g,e; r, J) 
a E {p)ai{q) 




(5.13) 



p + e q + e 



A direct comparison with equation (5.5) shows that for J = r 

R2(p,t;q,e;r,r) = R 2 (p,e;q,e;r). 

Computing the action of the volume and closing the trace with the h Sk operator we 
obtain 



Tr h Sh Vh Sh 




l>-> 

■j E E a e( r ) R 2{p,e;q,e;r,J) 

€,e,e— ±1 J—r,r-\-2e 
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OP 




© © 



\/w(p + e, q + e, r + e, 1) x 



Cr-n) 




We define 

A®(p,e;q,e;r) = ^ ^ y/w(p + e, q + e, r + e, 1) R 2 (p, e; q, e; r, r) 



e=±l 



CP 



r+e. 



x a e (r)- 




Cr±<) 



= X] 7 A/w(p + e,g + e,r + e, 1) #i(r,e) R 2 (p,e;q,e,r). 



e=±l 



We have finally obtained the action of H^ <2> on non planar 3-vertices 



Hi® |«(p,«,r)> = E E 



e=±l e=±l 



,4 <2) (p, e; g, e; r) 



+A®{q 7 e; r, e;p) 




+ A (2) {r,e;p,e; q) 



(5.14) 




(5.15) 



where the coeficients A 2 are given by (5.14). Notice that this is identical to the action 
of equation (5.11) with A® (p, e; q, e; r) replaced by A {2) (p, e; q, e; r). 



5.3. Other orderings 

Consider all the permutations of the four operator V, h Sk , /i~ and (h aij — h aji )/2. 
These are the "natural" orderings of the hamiltonian operator. Equations (4.19) and 
(4.20) exhaust the orderings in which the volume (the only operator that "grasps" 
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Table 1. Matrix elements A' 1 ' (p, e; q, e; r) and A® (p, e; q, e; r) of the non-symmetric 
and symmetric euclidean Thiemann's hamiltonian constraints, for the simplest 3- 
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new colors of tlic 


p 


q 


r 


e 
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A®(p,e;q,e;r) 




A< 2 )(p,e;g,f;r) 


1 





1 


1 


1 


(2,1,1) 





1/16^2 








1 


-1 


( NO ) 














-1 


1 


(0,1,1) 














-1 


-1 


{ J.NO ) 








1 


l 





1 


1 


(2,2,0) 
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-1 


( NO ) 
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( NO ) 
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-1 


( NO ) 














-1 


1 


(1,1,2) 
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-1 
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( NO ) 


o 





2 


2 





1 


1 


(3,3,0) 














1 


-1 


(3,1,0) 














-1 


I 




o 
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-1 


-1 


(1,1,0) 


o 


o 


2 


l 


1 


1 


1 
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1/16^2 
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1 


1 

- 1 
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_1 


I 


(1 2 1) 


5 /48 \/2 


^ /48 \/2 








-1 
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(1,0,1) 


1/32 v 7 ^ 





1 


2 


1 


1 


1 


(2,3,1) 


-1/16^2 


-1/16^5 








1 


-1 


(2,1,1) 


-5/48^2 


-5/48^2 








-1 


1 


(0,3,1) 














-1 


-1 


(0,1,1) 


-1/32^2 





1 


1 


2 


1 


1 


(2,2,2) 














1 


-1 


(2,0,2) 


1/24^12 











-1 


1 


(0,2,2) 


-1/24^12 











-1 


-1 


(0,0,2) 








2 


2 


2 


1 


1 


(3,3,2) 














1 


-1 


(3,1,2) 


-1/12(^5- ^2) 


-1/24-^30 








-1 


1 


(1,3,2) 


1/12(^5- ^2) 


1/24 ^0 








-1 


-1 


(1,1,2) 









the cylindrical function) sits between h Sk and h s . But these two are the only non 
vanishing cases, for the following reason. The h operators commute with each other 
and therefore we have to consider only the permutations of the three operator: V, 
O(o) = h Sk h~ k l and O^) — (h aij —h aji )/2. By the retracting identity, the O(o) operator 
is a two-index operator in the spin representation associated to the vertex v while 
the 0( 2 ) operator is a two-index operator in the spin 1 representation associated to 
the same vertex v. But in the hamiltonian constraint there is a trace, and thus in all 
the orderings of the operators we have to take a trace in the representation associated 
to the vertex. The volume operator modifies only the recoupling in the vertex and 
not the representation associated to it. Therefore we will always have a trace between 
a spin and a spin 1 representation, which gives zero. We can conclude that (4.17) 
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and (4.18) are the two only natural orderings for the HCO and their action reduces 
to the computation of equations (4.19) and (4.20). 

Following Thiemann p2j, we consider the symmetric euclidean HCO 



nf ymm [N\ = \ (h^[n} + h^[n?) . 

Using the transformation properties of h under conjugation, we have 

nf ymm [N] = \ (n^[N] + h^[n}) . 

The coefficients A s for the symmetric operators on trivalent vertices are 

A m (p, e;q,e; r)+A (2) (p, e; q, e; r) 
In table [l] we give the values of A (1) and A m for the simplest 3-vertices. 



A s (p,e;q,e;r) = - 



(5.16) 



(5.17) 



(5.18) 



6. Conclusions 

We have computed matrix elements of Thiemann's euclidean hamiltonian constraint, 
for its two natural orderings. In particular we have computed the action of the operator 
on any trivalent vertex in closed algebraic form. 

We collect here our results and express them in compact form. Following (45) , we 
define the operator D v jk& that acts on the vertex v by adding two new vertices on 
the edges J and K adjacent to v, and altering the colors of the two newly created 
edges by e, e = ±1 



D 



vJKle 





(6.1) 



Using this operator, the action of Thiemann's HCO over a generic spin-network state 
is 



H E [N]* s = il Q J2 N M 



1 



tiev(s) 



I,J,K£S(v) 
e,e=±l 



E(v) 



AvIJKU D v jKee 



(6.2) 



where V(s) is the set of the vertices of s and £(v) is the set of the edges adjacent to 
v. In general, A v jjK S g are matrices in the finite dimensional space of the contractors 
at the vertex i>, functions of the colors of the adjacent edges. On trivalent vertices, 
the space of these contractors has dimension one, and therefore the matrices A v ukzi 
are numbers. Let s have trivalent gauge-invariant vertices only. Denote as (r, p, q) the 
colors associated to the edges /, J, if of a vertex v. Then 



■ =^ E EE [A»(p,e;g,e;r)A,JX« 

veV(s) e=±l e=±l 

+A® (q, e; r, e; p)D vKIli + A® (r, e; p, e; q)D vIjrt 



(6.3) 
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where i — 1,2 for the two orderings we have studied. The explicit values of the matrix 
elements A (1) (p, e; g, e; r) and A (2> (p, e; g, e; r) are 



A m (p,e;q,e;r) = ^ - ^w(p, g, r + e, 1) i?i(r,e) R 2 {p,e;q,e,r). (6.4) 

e=±l 

A (2) (p,e;g,g;r) = ^ - + e, 9 + e, r + e, 1) i?i(r,e) R 2 (p,e;q,e,r). (6.5) 



where i?i(r, e), i?2 (p, e; g, e, ?") and w(a, 6, c, 1) are explicitly given by (5.8) (5.9) and 
(A. 17). We give some values of matrix elements in Table 1. 

The technique we have used can be applied in a straightforward manner to vertices 
of any valence. The only limitation is given by the lack of an explicit algebraic formula 



corresponding to equations ( A. 11 ) and (A. 17), for vertices of valence greater than four. 
The action of the volume operator is completly known; however, explicit values of its 
matrix elements are so far known throughout an algorithmic procedure Q only. It 
is interesting to note that nothing in the action of the operator seems to depend 
on the continuous moduli studied in [Q, which distinguish knots with high valence 
intersections. This might indicate that these parameters are superfluous and the 
theory lives actually in a smaller space. In this regard, see p7|. Finally, we expect 



that the technique used here can be applied to the second term in (4.7) as well, and 
therefore we expect that the total lorentzian hamiltonian constraint can be analyzed 
along the lines of this work. 
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Appendix A. Explicit formulas for the volume operator. 



A discussed in section 2.4, there exist two different versions of the volume operator 
[E5I. They have the same action on generic trivalent vertices. The volume operator 



acts on a cylindrical function \& 7 ,,/, 



13L \U 



by: 




(A.l) 




16-3! 



Y e(ei,ej,e K )Wij K 



ejue,;ue K =ii 



where, V("f) is the set of all vertices of 7, gi = ji(A(ei)), and Xi = X(gi) is the right 
invariant vector field on SU(2). All edges are chosen as outgoing. Notice that the 
action of the operator is given as a sum over terms in which three edges are "grasped" 
namely acted upon by one of the left-invariant vector fields. 
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'123 




Figure Al. Graphical representation of the grasping of the volume operator in the 
case of a gauge non-invariant 3-valent vertex. 



Here e(e/,ej,e/<) = sgn(det(ej(0), ej(0), ej<(0))) and we have denned 



W UK = 2 e„ fc XIX J T X 



(A.2) 



(q in jD| Q is related to W in || by W [UK] = 2 qu K .) As noted in Q, the 
difference between the volume introduced by Rovelli and Smolin in JTT| | and the one 
later introduced by Ashtekar and Lewandowski p3| is given by 



^ s _ 3! 



ejuejue K =D 
i 1 



16 W 



[7 J-ftT] 



16 3! e(e J ,e, 7 ,e K )W r [/ jK] 



ejuejue K =» 



(A.3) 
(A.4) 



This difference is due to a different regularization procedure used to derive the operator 
from its classical expression. From equations (A.4) and (A.3) it follows that the action 
of the two operators coincides when evaluated on generic (non coplanar) trivalent 
vertices, since there is only one term in the sum. 

Notice that the action of this operator is defined on gauge non-invariant states as 
well, since it is defined only in terms of left invariant vector field and therefore over 
the whole A space. Furthermore, explicit inspection of the definition of the operator 
shows that its action on a gauge non-invariant n-valent vertex coincides with one of 
the terms of its action on an a (n-fl)-valent vertex. 

To see this, consider the action of the volume on a 3-valent gauge non-invariant 
vertex in which the three colors of the adjacent edges are p,q,r, the color of the 
representation at the vertex is t, and the intertwiner is c. And consider the action 
of the volume on a gauge-invariant 4-valent vertex with edges colored p, q, r, t and 
intertwiner c. The latter is given by the sum of of four terms. In one of these four 
terms, the three left invariant vector fields in (A.2) grasp the three edges colored p, q, r. 
It is obvious that the resulting algebra is precisely the same as the one of the single 
term in the action of the volume on the trivalent gauge non-invariant vertex. 

Using this, we have, in graphical form, the key result on which the computations 
of this paper are based:: 



[012] 



J 



(A.5) 
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where we have indicated by [012] the triple of edges colored q,p and r. In the 
computations in the text, we have the particular case in which r is r ± 1 and t = 1. 
We write 



W m] I vj) = X) W$L/(q,p,r± 1,1) 



(A.6) 



where | Ur) is a compact notation for the vertex with internal edge colored /. Wfo^j 1 



was explicitly computed in in terms of 9-j and 6-j symbols: 

q p I 

qp (r± 1) A/ <( g p J V Tef 



it; 



(4) / 
[012] J 



(q,p,r± 1,1) 



qp (r±l) 



2 2 2 



J 7 1 
r ± 1 r± 1 2 



6(q,p,I)6(r± 1,1,1)0(2, J, I) 



r±L 



© © 







r ± 1 


2 







(A.7) 



This is the action on non-normalized states. Inserting the vertex normalization 
(equation (7.28) of ||) we get: 



r±l 



© © 



q _r±l 

© © 



[012; 



^ qp l«J>JV = 53iV^ [012] j / (p,g,r±l,l)^ 



Thus we can write the action of the (square of the) volume operator on the normalized 
states as 



W [m2] I vj) N = 53 {W [m2] /G>' ^ r ± l ) I ^)^' 



(A.8) 



where Wpj^ir'CP; 9: r ± 1, 1) are the matrix elements of the T4^ i2] vertex operator in 
the normalized basis. Their values are explicitly given by: 



W, 



(4) / 
[012] 



(p,q,r± 1,1) 



qp (r ±1) 



(A.9) 
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In this form it is obvious that W^qj^i j 1 (p, q,r ± 1, 1) is symmetric in I and J. 

An important observation is that I and J can take only two values in order to be 
compatible with r ± 1 and 1 at the triple intersection; these are r and r + 2_for the 
(+) case; and r and r — 2 in the (— ) case. But since the diagonal entries of Wj 1 are 
zero, there is only one possibility, which is I = r ± 2. Thus the eigenvalues of the 2x2 
matrix VT are ±i w(q,p, r ± 1, 1) with 



ttf(p, g, r ± 1, 1) = A / W$ 2]r r±2 (p, q,r±l, 1) 



(A.10) 
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From the eigenvalues of the W matrices we immediately obtain the eigenvalues of the 
volume operator. In fact, since the W matrix is a 2 x 2 antisymmetric matrix, its 
absolute value is a multiple of the identity. Thus the volume acts diagonally in the 
case in which we are interested, much simplifying the calculations. Taking the square 
root of the absolute value yields 



V v 



J 



in 



= -±y/w{q,p,r±l,l) 




r±i 



(A.ll) 



Let us compute the values of w. Starting from equation (A. 9) the matrix elements 

-(4) 
[012]- 



of wjk^J 1 are given by the chromatic evaluation: 



W/ [oi2]( a ' & ' c ' rf )"' / = a bc 



C£) C£ 



\ © © © © 
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b /* 


c 


2 


J 





(A.12) 



where the matrices W[q 12 i (a, b, c, d)j I are real antisymmetric matrices and have 
elements different from zero only if | J — 1\ — 2. Using equation (C.9) and applying 
(C.ll) to the two vertices v(a, 2, a), v(a, b, I) we obtain: 



W$ 2] (a,b,c, d)/ 



cp 



u y, i. L, 

\ ffi ffi S S 



(A.13) 
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and (C.14). Defining t = (I-J)/2 



This expression is evaluated using 

and e = (I — J)/2 we note that (a, 6, c, c?) j 7 is different from zero only if e = ±1 

and all the 3- vertices in equation (A.13) are admissible. In this way, we obtain the 
following explicit formula for the matrix elements (e = ±1): 

1 a+b+t+3c+d +t+3 



W$ 2] (a,b,c,d) t +l= -e(-l)- 
1 + a + b 



_4i(i + 2) 2 
tl + a + t-bl + b + t-a 



1+c+d-tl 



t-dl 



2 

d + t-c 



(A.14) 



2 2 2 

This formula can be used to obtain explicit expressions for the eigenvalues of the 

' (4) 
[012 



iW,n}m{a,b,c,d)\ operator when they are 2x2 matrices, as in the case of equation 



(A.ll 



w{a, a, 1, 1) = 



a(a 



(A.15) 
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„ cl b c (a + b + c) , . . 

w(a,b,c,a + b + c-2) = 2W ^ (A.16) 

1 l {a + b + c + 3){l + a + b-c)(l+c-a + b)(l + c + a-b) 
w(a,b,c,l) = -\J — (A.17) 

Appendix B. Graphical calculus in the connection representation. 

Let us recall the basic elements of Penrose binor calculus , which is the basis for 
expressing the cylindrical functions of generalized connections in graphical form |37j . 
The key idea of Penrose is to rewrite any tensor expression in which there are sums 
of dummy indices in a graphical way [p4[ . Penrose represents the basic elements of 
spinor calculus (i.e., tensor expression with indices A, B, ... — 1,2) as 



°c 



VA 



(B.l) 



TaSb 



and assigns a minus sign to each crossingf , i.e: 



x£B = 6f6 B = - ^ (B.2) 

Using this rule it is possible to represent any SU (2) tensor expression in a graphical 
way as follows: (1) Define the "up" direction in the plane; (2) Draw boxes with the 
name of the SU(2) tensor (except for the S and the e that are simply represented 
by lines) with as many slots going up as the number of contravariant indices and as 
many slots going down as the number of covariant indices; (3) For each dummy index 
connect with a line the corresponding slots of the boxes; (4) Assign a "i" factor to 
each minimum or maximum of the lines; (5) Assign a minus sign to each crossings of 
lines. Conversely, any curve can expressed as a product of (Ts and e's. 

An unexpected feature of this graphical calculus is its way of expressing the 
operation of taking the trace. Since any contravariant index of a tensor is represented 
more "up" than any covariant index of the same tensor, there is no way to write 
directly the operation of taking the trace unless one uses the antisymmetric 2-tensors 
6 A B and e AB as TrXg = Xg5 B = eA' At B ' B X^S^, . This becomes, graphically, 



TrX A B = S B X A B = - \g} . (B.3) 
Moreover, the basic binor identity is graphically given by 

f I | + X = $b 5a + 8f + (-1) e AB e CD = 0, (B.4) 



perhaps the most relevant equation of loop quantum gravity. 

In particular, the irreducible representation of color n can be constructed as the 
symmetrization II„ (in Penrose's graphical representation an anti-symmetrization 



because of equation (B.2)) of the tensor product of n fundamental representations 
^ B — U(g) B {g being the group element). Denoting by 11^ the normalized 



f The binor representation is a graphical representation in the plane. 
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symmetrizer (graphically an anti-symmetrizer 
representation is given by: 



nWp n = l^(-i)W 



see p7j]), the color n irreducible 



(B.5) 



J/(«/) = 



(B.6) 



(P„ represents rt parallel lines, P„ is the graphical representation in terms of lines of 
the permutation p, \p\ is the parity of the permutation and a line labeled by a positive 
integer n represents n non- intersecting parallel line). 

Analogously, it is possible to construct an explicit graphical representation of the 
Clebsch-Gordan intertwining matrix i.e., of the matrix that represents the coupling of 
3 (or more) irreducible representations of the SU(2) group. The representation of the 
Clebsch-Gordan intertwining matrix in the binor formalism is given by the special sum 
of "tangles" denoted as the 3-vertex. Each line of the vertex is labeled by a positive 
integer a, b or c and is defined as: 



4ef 




(a + b-c)/2 
(b + c-a)/2 
(c + a-b)/2 



(B.7) 



where m,n,p are positive integers. This last condition is called the admissibility 
condition for the 3-vertex. By the Wigner-Eckart theorem, or more precisely, by 
the version of it due to Yutsin, Levinson and Vanagas any invariant tensorial 
intertwining matrix representing the coupling of n representations of a compact group 
can be expressed as the product of Clebsch-Gordan coefficients or, which is the same, 
in terms of a trivalent decomposition. It is easy to see that any contractor c a of a 
given vertex v a with n incoming (outgoing) edges of color Pq,. . . ,P„_i, can be written 
in terms of a linear combination of the trivalent contractors obtained ordering the 
edges and assigning n-3 compatible integers, i.e. by a trivalent decomposition of the 
vertex. Graphically: 



P? 



Pi— a 



Pi 
P, 



U Pn 

, Pn 



Pi-3 

Pi -Y* 2 *3 i„_K_ P n -2 ' ( R8 ) 
P / \ Pn-1 

n _2 that satisfy the admissibility 




where the sum is extended to all the integers 12, ■ • • 
conditions for the 3- vertices. 

An important characteristic of the binor representation, related to the existence of 
the natural "metric" (an invariant way to relate contravariant and covariant indices) 
eab in the 2-spinor space, it is its independence from the orientation chosen for 
the edges. In fact, the SU{2) relation {/(.g -1 )^ = e-Ac eBD U{g) c D is graphically 
represented (in the binor representation) by: 



-14 



while 



(B.9) 
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is the related identity in case of the n-representation. 

In this graphical notation the left invariant vector fields associate to an edge are 
represented by: 




(B.10) 



From the identity (that can be proven in a straightforward way): 

e ufe ^® =2 

we have the important result 

in/ mi irijc 

|*^|«*=^ J" I | (B.12) 



2 




(B.H) 



which is the key "bridge" showing that the volu me op erator that was defined in the 
loop representation in [O, based on the r.h.s. of ( B.12|) , is indeed the same operation 
as the volume that was defined in the connection representation in [fl3|| , based on the 
l.h.s. of JB.12(). 



Appendix C. Tangle-theoretic recoupling theory. 

One of the main results of the recoupling theory of colored knots and links with 
trivalent vertices |u| is the computation of the Kauffman bracket for framed spin 
networks. The framing refers to the fact that in the computations one keeps track of 
over- and under-crossing. In our context this distinction has no meaning. Therefore 
we use a simplified version of the recoupling theory, namely we replace the deformation 
parameter A by its "classical" value —1, relevant to our case. This case corresponds 
to the standard recoupling theory of the SU(2) group. In the formulae from |31f| , we 
thus replace the deformed (or quantum) integers by ordinary integers. We list here 
the basic formulae which we use. For more details see |3l| . This beauty of the theory 
is that it is completely determined by the following two tangle properties: 

(X ) = MX ) + I | >, (c.i) 

( O U K ) =d ( K ), (C.2) 

where d — —A 2 — A~ 2 and K is any diagram that does not intersect the added loop. 
The graphical representation of the states that we have defined in section 3 does satisfy 
these two equations (with A = — 1) and therefore all recoupling theory equations can 
be applied. 

There is a subtle point that needs to be clarified here. When tangle theoretical 
recoupling theory is defined on a pane, it corresponds, via Penrose binor calculus to the 
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algebra of e^ B and tensors only. The states in the connection representation are of 
course not given just by these tensors, because they contain the parallel propagators. 
In Penrose binor calculus, these tensors are represented as boxes. Thus, tangle 
theoretical identities can be applied only in between these boxes - for instance, a 
closed line cannot be replaced by a factor d, as in (C.2), if it includes a box. However, 
in the graphical representation of section 3, we have coded the information of the 
presence of the parallel propagators in the fact that a line runs along an edge. A short 
reflection will convince the reader that if one applies recoupling theory to the tangles 
defined on the extended graph (which is not a plane, but rather an oriented 2-d surface 
with a complicated topology) that is, one defines tangles as equivalence classes under 
smooth deformations on the extended graph, then this is equivalent to respecting the 
presence of the boxes - for instance, a generic line will not be contractible to a circle, 
and therefore equation (C.2) can be applied only within the vertices, and not for a loop 
wrapped around the extended graph. This is why recoupling theory can be rigorously 
applied to tangles defined on the extended graphs. 

One of the consequences of equation (C.2) is that any closed loop that can be shrink 
to a point within the extended graph can be replaced by its loop value, which is equal 
to —2. Indeed, any closed tangle that shrinks to a point reduces to a number and the 
operation of computing this number is usually referred as the chromatic evaluation. 
The most useful evaluations are the following: 
(1) The symmetrizer 



A„ = 




= (-!)> +1). 



(C.3) 



(2) The 9 net 



6(a, b, c) = 




(-l) m+n +P(m + n + p + 1)! to! n\ p\ 
al bl c! 



(C.4) 



where m = (a + b — c)/2, n = (b + c — a)/2, p = (c + a — b)/2. 
(3) The Tetrahedral net 



Tet 



ABE 
CDF 



where 



Oi 

a 4 




A + D + E 

2 ' 
B + C + E 



A + B + F 



C + D + F 



(-l) s (S + l)\ 



, U 1I, (s-aty. n ; (bj-sy. 

B + D + E + F 



(C.5) 



bi 
b 2 



A+C+E+F 



A+B+C+D 



to = max{ai}, M — mm{bj}, 

£ = A\B\C\D\E\F\, 1 = II,,;',, -/,:■! 
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The most important result of the theory is the recoupling theorem which states: 

b c b \y c 

,/C =E{° J j} X («) 

d 1 ad 



where the quantities are the s«(2) 6 — j symbols (normalized as in and their 
explicit values, in terms of chromatic evaluations, are: 

Ai Tet 



a b i 
c d j 



a b i 
c d j 



6(a, d, i)0(b, c, i) 



(C.7) 



From the previous formulae and definitions other relations can be easily derived. 
In particular, exchanging the ordering of the three line in a 3-vertex yields a sign 
factor: 




a b 



= XT 



(C.8) 



where Xf = (-l)(°+ 6 - c )/ 2 (_l)(«'+6'-c')/2 j and x > = x ( x+ 2); The following reduction 
holds 



fc. 
a . 

i 



c\d 




(C.9) 



Another key formula that can be derived is the one we use for evaluating the action 
of the parallel propagator operator on spin network states: 



nil Ui a Kl 1 

£o.(n) n+e = „ + i nr-i 

E =±l |n|l A ' 



ffc — 

nT| ] 



(CIO) 



Also 



P 9 
2 



P x 2 ? 



(C.ll) 



Finally, using equations (C.6) and (C.7), we have the chromatic evaluations (e = 0, ±1) 
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(C.13) 




t±^ _ (3 + a + b + t)(l + o + 1 - b)(l + b + t - a) 




4t(i + l)(l + a + &-i) 



These formulas are sufficient for the computations in the paper. 
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